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Anomalously large oxygen-ordering contribution to the thermal expansion of untwinned 
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We present high-resolution capacitance dilatometry studies from 5 - 500 K of untwinned YBa2Cu30j^ 

(Y123) single crystals for x ~ 6.95 and x = 7.0. Large contributions to the thermal expansivities due to 
0-ordering are found for x ~ 6.95, which disappear below a kinetic glass-Uke transition near room 
temperature. The kinetics at this glass transition is governed by an energy barrier of 0.98 ± 0.07 eV, in very 
good agreement with other 0-ordering studies. Using thermodynamic arguments, we show that 0-ordering 
in the Y123 system is particularly sensitive to uniaxial pressure (stress) along the chain axis and that the 
lack of well-ordered chains in Ndl23 and La 123 is most likely a consequence of a chemical-pressure effect. 

PACS number: 74.72.Bk, 64.70.Pf, 74.62.Fj, 65.50.-i-m 



The superconducting and normal state properties of 
YBa2Cu30jj^ (Y123) depend strongly both on the total 

oxygen content (O-content) x and on the degree of O-order 
in the CuO chain layer, and a great deal of experimental and 
theoretical effort has gone into studying both the different 
types of O-ordered structures and their correlation with 
superconducting properties [1-23]. Most of these studies 
have concentrated on samples with O-contents in the range 
6.3 < X < 6.7, probably because superconductivity is 
particularly sensitive to the degree of O-order in this range. 
The kinetics of O-rearrangements occurs via a thermally 
activated diffusion process with a time constant on the order 
of hours at room temperature, and many interesting effects 
have been observed by room-temperature annealing studies 
[7-13]. 

There are however only a few studies in which 0-ordering 
effects have been reported near the optimal doping level 
(6.9<x<7.0), where there are only a few O-vacancies and 
much smaller ordering effects are expected. Nevertheless, 
resistivity studies of Lavrov et al. clearly demonstrated that 
O-ordering occurs even in nearly fully doped (x~6.95) 
samples [14]. Further evidence for O-ordering in optimally 
doped crystals has also been obtained by high-temperature 
quench experiments of Erb et al., in which it was clearly 
shown that O-vacancy clusters form in slowly cooled 
crystals and that these clusters form strong flux-pinning 
centers, which are responsible for the 'fish-tail' effect [16]. 

In this Letter we present a thermodynamic and kinetic 
study of O-ordering processes in near optimally doped 
(x~6.95) Y123 untwinned single crystals using high- 
resolution thermal expansion measurements. Surprisingly 
large contributions to the linear thermal expansion 
coefficients (expansivities) from O-ordering are observed. 
These are, as expected, absent in the fully oxygenated 
(x=7.0) state. The large anomalies imply that O-ordering is 
highly sensitive to slight dimensional changes of the crystal 



lattice. The largest contribution occurs along the chain axis 
(b-axis), from which the peculiar result, that compressing 
the chain axis increases the degree of O-order, can be 
inferred. Our results clearly show that the lack of well- 
ordered chains in Ndl23 and Lal23 [17] is a consequence 
mainly of the larger (in comparison to Y123) unit cell. 

Our investigations have been performed on Y123 single 
crystals grown by a pulling technique, detwinned at a 
uniaxial pressure of 10 MPa and oxygenated at 490°C [24]. 
This procedure results in a slightly overdoped sample with 
an O-content x of approximately 6.95. In order to get to the 
fully oxygenated state with x = 7.0, the crystals were 
annealed at 376 bar O-pressure and 400°C for 140 hours. 
The same crystal, with linear dimensions of 6 x 3 x 2 mm^ 
for the different orthorhombic axes a, b and c, respectively, 
as used in a previous thermal expansion investigation [25] 
was used presently. The thermal expansion was measured 
with two different high-resolution capacitance dilatometers 
in the range 5 - 300 K and 150 - 500 K respectively, using 
constant heating and cooling rates between 2.5 and 20 mK/s. 
Data points were taken every 100 mK. 
In Fig.l we present the linear thermal expansivities, a(T) = 

X*"^, of Y 123 along the three orthorhombic axes in the 

temperature range 5 - 500 K for two different x values [x ~ 
6.95 (solid line) and x = 7.0 (dashed line)]. In addition to 
the well-studied superconducting anomalies at T^ [25-27], a 
quite large and broad a(T) anomaly is observed for the 
X ~ 6.95 sample near room temperature, or more precisely at 
Tg = 280 K, along each of the axes. No comparable effect is 
seen for x = 7.0, and in the following we show that this 
anomaly results from a kinetic glass-like transition due to 
O-ordering. 

The shape of this anomaly (width [28] and hysteresis for 
heating and cooling curves) are typical of kinetic 'glass' 
transitions [29-31]. A comparison between the x ~ 6.95 and 
X = 7.0 data sets clearly indicates that, as expected [29-31], 




Fig.l. Expansivity versus temperature of YBa2Cu30x for the a, b 
and c axes for x ~ 6.95 (solid lines) and x = 7.0 (dashed 
lines). In addition to the superconducting phase transition 
at Tc, a glass-like transition can be observed at Tg in the 
data with x ~ 6.95. The arrows indicate heating and cooUng 
cycles. 

there are additional contributions (Aa) to the expansivities 
above Tg, which are absent below Tg. The contribution 
above Tg is positive for the b and c-axis, and negative for 
the a-axis [32] (see Fig. 1). The largest effect is found for 
the b-axis, for which we have also examined the expansivity 
as a fiinction of cooling rate for -2.5, -5, -10 and -20 mK/s. 
As seen in Fig. 2, the glass transition shifts to higher 
temperatures with increasing cooling rates, which clearly 
indicates the kinetic nature of the transition. Here we have 
subtracted a background, which was constructed by shifting 
the data of the b-axis for x = 7 by a small value, from the 
original data. In order to extract a well-defined Tg from 
these broad transitions, we have normalized the data in 
Fig. 2 by the equilibrium expansivity Aab"''"'' (dashed line). 
Below Tg, Aab"*"' is obtained by an extrapolation of the 
high-temperature data to lower temperatures. In Fig. 3 these 
normalized curves are plotted versus -T"' in order to show 
that the glass transition shifts in temperature, while the 
shape of the glass-transition does not change. Tg values can 
now easily be extracted from Fig. 3 as the midpoints of the 
transitions, and these are shown in an Arrhenius plot (see 
inset in Fig. 3). The data show a good linear dependence, 
indicative for a thermally activated relaxation time [29-31] 



T(r) = v„-'-e 



k„T 



with an activation energy Ea 
relaxation time x at T„ 



(1) 

0.98 ± 0.07 eV. The 



(2) 



can be directly related to the coohng rate, r, in our 
experiment [33], yielding the pre-factor Vq = I-IO'^ Hz in 
Eq. 1. Our value for the activation energy E^, which is 
physically related to the barrier for O-hopping, is in very 
good agreement with those found by other kinetic studies of 
0-ordering [5,8,15,34]. This and the fact, that the glass 
transition vanishes for x = 7, provides clear evidence, that 
the expansivity anomahes in the x ~ 6.95 sample result from 
the freezing of O-ordering processes. This ordering occurs 
on a 'local' scale, as indicated by the short extrapolated 
diffusion lengths using the experimental diffusion constants 
[34] and our relaxation time t(280 K) = 350 s. We stress 
that we see no indications for the so-called '240 K - phase 
transition', which has also been discussed in connection 
with O-ordering [35]. In the following, we analyze the 
thermodynamics of the O-ordering above Tg, where thermal 
equiUbrium is guaranteed. 

The additional contribution to the expansivity Aa above 
Tg (see Fig. 2) is a direct meastire of the uniaxial pressure 
(stress) dependence of the degree of O-order at constant 
temperature [36] 

1 dS, 

Aa.. = 



^ Oxygen 



(i = a, b, c). 



(3) 



Here, Soxygen is the molar entropy associated with O-disorder 
and Vmoi is the molar volume. Inserting the measured Aa 
values at 340 K (Aaa = -5-10"^K-', Aat = 30-10"^K"\ 
Aa^ = 8-10"^K"') into Eq. 3, we obtain dSoxygen/dpa = 
56 n]J/(mol-K-GPa), dSoxygen/dpb = -311 mJ/(mol-K-GPa) 
and dSoxygen/dpc = -85 mJ/(mol-K-GPa). This shows that 
uniaxial presstire (stress) along the a-axis decreases the 
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Fig.3. Normalized plot of Aa(,(T) versus -T for x ~ 6.95. The 
inset shows an Arrhenius-plot of the cooling rate and 
relaxation time versus Tg"' (see text for details). 

O-order, whereas uniaxial pressure along the b- and c-axis 
increases the O-order. Hydrostatic pressure is expected to 
increase the O-order since dSoxygen/dphyd = S dSoxygen/dpi = 
-0.34 J/(mol-K-GPa). We can also calculate the change of 
Soxygen resulting from volume changes 

AS 0-34 — fr- 

=B-AP-V^,1 = „ (4) 



volume change 

using the bulk modulus B = 100 GPa [37] and the 
O-ordering contribution to the volume expansivity Ap. In 
order to get a feeUng for the magnitude of these pressure 
(volume) effects, it is instructive to calculate the maximimi 
entropy difference ASg^gen'-*^'^^^'^ totally ordered and 
totally disordered states. For x = 6.95 (5 % 0-vacancies) 
and assuming that the disorder is limited to the disorder in 
the chains (no 0-5 site occupations), AS ^^^^ can be easily 

calculated using elementary statistical mechanics to equal 
1.65 J/(mol-K) [38]. Thus, a relatively small (1 %) volume 
increase (decrease) will result in a quite large (20 % of 
AS Oxygen ^ decrease (increase) of Soxygen- This result, which 

is a direct thermodynamic consequence of our data, provides 
a simple explanation for the different degrees of O-order 
observed in R 123 (R= rare-earth) compounds [17]. In R 123 
compounds with a larger (smaller) unit cell (in comparison 
to Y123) the degree of O-order should be much lower 
(much higher) than in Y123. Ndl23 e.g. has a unit-cell 
which is about 2.6 % larger than that of Y123 [39] and 
should therefore have an entropy, which is about 
0.88 J/mol-K larger than in Y123 - this entropy-difference is 



about 50 % of ASo^ggn , and suggests that Ndl23 is nearly 

fully disordered, in agreement with experiment [17]. Also, 
the failure to detect any glass-transition in Ndl23 with 
x ~ 6.95 up to 500 K using dilatometry [18], which implies 
that dS/dV = (see Eq. 4), suggests that Soxygcn is at its 
maximum value. Ybl23, on the other hand, has a smaller 
unit cell and is expected to have a higher degree of chain 
order. Our results suggest that the difference between the 
degree of O-ordering in the R-123 can be explained by a 
simple chemical-pressure effect. 

The additional contribution Attb (see fig.2) and, thus, 
dS/dpb decreases strongly in magnitude with increasing 
temperature, and this provides additional thermodynamic 
information about the O-ordering. An analogous 
contribution has to also occur in the specific heat Cp, but, 
unfortunately, even high-resolution measurements have not 
detected this anomaly [40]. These specific heat data, 
however, provide an upper limit of about 2 J/(mol K) for 
this contribution at 300 K. It is instructive to compare this 
value with the difference AC between the specific heats at 
constant volume Cy and constant pressure Cp, which can be 
directly calculated from our volume expansivity P using 

(5) 



AC = C -C =T V 



BP'. 



Just above Tg, AC = 0.8 J/(mol-K), which is about 40 % of 
the maximally estimated Cp contribution. This indicates 
that, in addition to the usual increase in entropy due to 
increasing temperature, the volume change due to the 
thermal expansion plays a major, if not dominating, role in 
the increase of entropy with increasing temperature; this is 
because the dS/dV term is so large. Another way to quantify 
the strong volume dependence is through a partial Griineisen 
parameter for O-ordering 

. .O-ordering ^mol ^ 

/ Griineisen AC ' 

Inserting the maximum estimate for the specific heat 
anomaly just above Tg into Eq. 6 yields a lower limit for 
Yoruneisen (> 17). Au uppcr limit (YcrQnciscn < 43), ou the Other 
hand, is obtained by reaUzing that ACv must be greater than 
or equal to zero, from which a lower limit of ACp > 
0.8 J/(mol-K) can be inferred (see Eq. 5). Such large values 
of YGriineisen clcarly document the enormous pressure/volume 
dependence of the O-ordering effects in Y 123. 

So far we have not discussed either the spatial 
arrangements or the microscopic origin of these ordering 
processes. From optical measurements it is clear that the 
chain length is growing as one lowers the temperature [19]. 
This is also compatible with the hole transfer from chains to 
planes, which is expected as the chains increase in length 
[20], observed by Lavrov at al. [14], as well as with the 
formation of O-clusters [16]. Theoretically, O-ordering has 
been modeled using the asymmetric-next-nearest-neighbor- 
interaction (ASYNNNl) model, in which several 
phenomenological parameters (VI > V2 > V3 > ...) are 
introduced to describe the O-interactions in the basal plane 



[21-23]. We note that we do not observe the thermodynamic 
phase-separation line (ortho-I ortho-I + ortho-II) 
predicted by this model [22] in our experimental 
temperature window (280 K < T < 500 K). The parameter 
V2 « 2000 K determines the length of the chains in this 
model. A recent theoretical study has argued that a single 
V2 parameter is too simple, and that this interaction depends 
on chain length as well as on x [41]. Our results show that 
V2 depends also strongly on the size of all lattice 
parameters, which demonstrates the complexity of this 
interaction. We should point out that this cannot be 
understood in terms of a simple steric effect involving just 
the b-axis parameter, since V2 also depends on the size of 
the other axes, and any ab-initio calculations of V2 will 
most likely have to include the electronic structure of the 
whole unit cell, including effects such as charge transfer etc. 

In conclusion, our thermal expansion data of untwinned 
single crystals clearly reveal the kinetics and, for the first 



time, also the thermodynamics of O-ordering effects in 
YBa2Cu306.95. These ordering effects manifest themselves in 
terms of a kinetic glass-like transition neiir room 
temperature. The large size of the glass-transition anomalies 
along all axes are directly correlated with the extreme 
dependence of the degree of O-order with the lattice 
parameters, and this provides a natural explanation of the 
different O-order (e.g.) in Ndl23 and Lal23 in terms of a 
simple chemical-pressure effect. Studies of O-ordering in 
R123 compounds using high-pressure or strain induced by 
lattice mismatch in thin films should provide further 
valuable information [42]. Theoretical modeling of the O- 
ordering via lattice gas models such as the ASYNNNI 
model, should incorporate this new results for a more 
quantitative comparison with experiments. 
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